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Abstract 

An algebraic transformation of the DeTemple-Wang half-integer approximation 
to the harmonic series produces the general formula and error estimate for the Ra- 
manujan expansion for the nth harmonic number into negative powers of the nth 
triangular number. We also discuss the history of the Ramanujan expansion for the 
nth harmonic number as well as sharp estimates of its accuracy, with complete proofs, 
and we compare it with other approximative formulas. 
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1 Introduction 



1.1 The Harmonic Series 

In 1350, Nicholas Oresme proved that the celebrated Harmonic Series, 
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is divergent. He actually proved a more precise result. If the n th partial sum of the harmonic 
series, today called the n th harmonic number, is denoted by symbol H n : 
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then Oresme proved that the inequality 
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holds for k = 2, 3, ... . This inequality gives an estimate of the speed of divergence. 

Almost four hundred years passed until Leonhard Euler, in 1755 [3] applied the 
Euler-Maclaurin sum formula to find the famous standard Euler asymptotic expansion 
for H n , 
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where B\~ denotes the k th Bernoulli number and 7 := 0.57721 • • • is Euler's constant. This 
gives a complete answer to the speed of divergence of H n in powers of -. 

Since then many mathematicians have contributed other approximative formulas for H n 
and have studied their velocity of divergence. We will present a detailed study of such a 
formula stated by Ramanujan, with complete proofs, as well as of some related formulas. 

1.2 Ramanujan's Formula 

Entry 9 of Chapter 38 of B. Berndt's edition of Ramanujan's Notebooks [21 p. 521] reads, 
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-, where n is a positive integer. Then, as n approaches infinity, 
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We note that m := is the nth triangular number, so that Ramanujan's expan- 

sion of H n is into powers of the reciprocal of the n th triangular number. 

Berndt's proof simply verifies (as he himself explicitly notes) that Ramanujan's expansion 
coincides with the standard Euler expansion (jl.4p . 

However, Berndt does not give the general formula for the coefficient of -\ in Ramanu- 
jan's expansion, nor does he prove that it is an asymptotic series in the sense that the 
error in the value obtained by stopping at any particular stage in Ramanujan's series is less 
than the next term in the series. Indeed we have been unable to find any error analysis of 
Ramanujan's series. 

We will prove the following theorem. 



Theorem 1. For any integer p ^ 1 define 
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where B 2 k{x) is the Bernoulli polynomial of order 2k. Put 
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where n is a positive integer. Then, for every integer r ^ 1, there exists a O r; < © r < 1, 
for which the following equation is true: 
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We observe that the formula for R p can be written symbolically as follows: 

i (ab 2 -iy . . 

R ^-2 V \—)^ (L9) 

where we write B 2m (^) in place of B 2m after carrying out the above expansion. 

We will also trace the history of Ramanujan's expansion as well and discuss the relative 
accuracy of his approximation when compared to other approximative formulas proposed by 
mathematicians. 



1.3 History of Ramanujan's Formula 

In 1885, two years before Ramanujan was born, Cesaro [1] proved the following. 

Theorem 2. For every positive integer n ^ 1 there exists a number c n , < c n < 1, such 
that the following approximation is valid: 



H n = - ln(2m) + 7 + — 
2 v 1 ' 12m 
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This gives the first two terms of Ramanujan's expansion, with an error term. The method 
of proof, different from ours, does not lend itself to generalization. We believe Cesaro's paper 
to be the first appearance in the literature of Ramanujan's expansion. 

Then, in 1904, Lodge, in a very interesting paper [8J, which later mathematicians inex- 
plicably (in our opinion) ignored, proved a version of the following two results. 

Theorem 3. For every positive integer n, define the quantity X n by the following equation: 
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Theorem 4. For every positive integer n, define the quantity A n by the following equation: 
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where < S„ < — . The constants in the expansion of A n all are the best possible. □ 
4042500 

These two theorems appeared, in much less precise form and with no error estimates, in 
Lodge [Sj. Lodge gives some numerical examples of the error in the approximative equation 
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in Theorem [31 he also presents the first two terms of A n from Theorem HI An asymptotic 
error estimate for Theorem [3] (with the incorrect constant ^ instead of ) appears as 
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Exercise 19 on page 460 in Bromwich [3]. 

Theorem and Theorem [7] are immediate corollaries of Theorem [ZJ 

The next appearance of the expansion of H n , into powers of the reciprocal of the n th 

triangular number, m = n ^ n+1 ^ , is Ramanujan's own expansion (11.51) . 
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1.4 Sharp Error Estimates 



Mathematicians have continued to offer alternate approximative formulas to Euler's. We 
cite the following formulas, which appear in order of increasing accuracy. 
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Formula 1 is the original Euler approximation, and it overestimates the true value of H n 

by terms of order -. 

J 12n 2 t 

Formula 2 is the Toth-Mare approximation, see [9], and it underestimates the true 

value of H n by terms of order -. 

3 72n 3 

Formula 3 is the Ramanujan-Lodge approximation, and it overestimates the true 

19 

value of H n by terms of order — — -, see [TO] . 

3150[n(n + 1)J 3 

Formula 4 is the DeTemple-Wang approximation, and it overestimates the true value 
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of H n by terms of order =— ■, see IB]. 
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In 2003, Chao-Ping Chen and Feng Qi [5] gave a proof of the following sharp form 
of the Toth-Mare approximation. 



Theorem 5. For any natural number n ^ 1, the following inequality is valid: 

(1.12) 




The constants 2 = .3652721 • • • and d\ are the best possible, and equality holds only 

1 — 7 6 

for n = 1 . □ 

The first statement of this theorem had been announced ten years earlier by the editors 
of the "Problems" section of the American Mathematical Monthly, 99 (1992), p. 685, as part 
of a commentary on the solution of Problem E 3432, but they did not publish the proof. So, 
the first published proof is apparently that of Chen and Qi. 
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In this paper we will prove new and sharp forms of the Ramanujan-Lodge approximation 
and the DeTemple-Wang approximation. 

Theorem 6 (Ramanujan-Lodge). For any natural number n ^ 1, the following inequality 
is valid: 
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holds only for n — I. 

Theorem 7 (DeTemple-Wang). For any natural number n ^ 1, the following inequality is 
valid: 
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DeTemple and Wang never stated this approximation to H n explicitly. They gave the 
asymptotic expansion of H n , cited below in Proposition [U and we developed the correspond- 
ing approximative formulas given above. 

All three theorems are corollaries of the following stronger theorem. 

Theorem 8. For any natural number n ^ 1, define f n , \ n , and d n by 
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respectively. Then for any natural number n ^ 1 the sequence {f n } is monotonically 
decreasing while the sequences {A n } and {d n } are monotonically increasing. 

Chen and Qi [5] proved that the sequence {f n } decreases monotonically. In this paper 
we will use their techniques to prove the monotonicity of the sequences {A n } and {d n }. 
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2 Proof of the sharp error estimates 
2.1 A few Lemmas 

Our proof is based on inequalities satisfied by the digamma function 



:= -^lnr(x) = ^4 
K J dx K J T(x) 
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which is the generalization of H n to the real variable x since ty(x) and if n satisfy the equation 
[H (6.3.2), p. 258]: 

#(n + l) = # n - T . (2.2) 

Lemma 1. For every x > i/iere ercist numbers 9 X and Q x , with < 9 X < 1 and < Q x < 1, 
for which the following equations are true: 
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Proof. Both formulas are well known. See, for example, [3 pp. 124-125]. 

Lemma 2. The following inequalities are true for x > 0: 
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Proof. The inequalities (12.51) are an immediate consequence of (I2.3P and the Taylor expansion 
of 
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which is an alternating series with the property that its sum is bracketed by two consecutive 
partial sums. 

For (12.61) we start with (12.41) . We conclude that 
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Now we multiply to all three components of the inequality by 2 and add 
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Lemma 3. The following inequalities are true for x > 0: 
111111 1 
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Proof. Similar to the proof of Lemma [2j □ 

2.2 Proof for the Ramanujan— Lodge approximation 

Proof of Theorem^ for {A n }. We solve (11.151) for A n and use ( 12.21) to obtain 

A n = : — 6n(n + 1). 

+ 1) - In y/n(n + 1) 
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A x := - - — 3x(x + 1), 

2#(x + l) -lnx(x + l) 

for all x > 0. Observe that 2A n = A n . 

We will show that that the derivative A' x > for x > 28. Computing the derivative we 
obtain 
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(by the remainder theorem), which is obviously positive for x > 28. Thus, the sequence 
{A n }, n ^ 29, is strictly increasing. Therefore, so is the sequence {A n }. 
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For n — 1, 2, 3, . . . , 28, we compute A n directly: 



1.1831718 
1.1961868 
1.1983668 
1.1990988 
1.1994300 
1.1996073 
1.1997131 

Therefore, the sequence {A n }, n ^ 1, is a strictly increasing sequence. 
Moreover, in Theorem [31 we proved that 
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2.3 Proof for the DeTemple— Wang Approximation 

Proof of Theorem^ for {d n }. Following the idea in the proof of the Lodge-Ramanujan ap- 
proximation, we solve (11.161) for d n and define the corresponding real-variable version. Let 
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We compute the derivative, ask when is it positive, clear the denominator and observe that 
we have to solve the inequality: 

{^T + !)} -48(x+i){^(x + l)-ln(x + i)} 2 > 0. 

By Lemma El the left hand side of this inequality is 
11111 1 



x + 



48 



1 1 

x 2x 2 



1 1 

6x 3 30x 5 



+ 



42x 7 
23 



11 



+ 



143 



^24x 2 2Ax 3 960x 4 160x 5 8064x 6 896x 7 30720x 8 y 
for all x > 0. This last quantity is equal to 
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The denominator is evidently positive for x > and the numerator can be written in the 
form 

p(x)(x — 4) + r, 

where 

p(x) = 548963242092 + 137248747452x + 34315688832a; 2 + 8564093760x 3 + 2138159872x 4 
+ 566849792x 5 + 111820800x 6 + 11547648x 7 + 178139136x 8 + 178139136x 9 , 

with remainder r = 2195843950359. 

Therefore, the numerator is clearly positive for x > 4, and therefore, the derivative d' x is 
also positive for x > 4. Finally, 

di = 3.73929752 • • ■ 
d 2 = 4.08925414 • • ■ 
d 3 = 4.13081174- •• 
di = 4.15288035- ■• 

Therefore {d n } is an increasing sequence for n ^ 1. 

Now, if we expand the formula for d n into an asymptotic series in powers of f , we 

obtain 
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(this is an immediate consequence of Proposition [1] below) and we conclude that 
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3 Proof of the general Ramanuj an— Lodge expansion 

Proof of Theorem [3 Our proof is founded on the half- integer approximation to H n due to 
DeTemple and Wang [6]: 

Proposition 1. For any positive integer r there exists a 6 T , with < 9 r < 1, for which the 
following equation is true: 



H n=H n +h)+1 + ^ I ^ + e r . I - T ^_, (3,) 
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and where B 2p (x) is the Bernoulli polynomial of order 2p 
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Since (n + |) 2 = 2m + |, we obtain 
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Substituting the right hand side of the last equation into the right hand side of (13. ip we 
obtain 

F n = ln(n + \) + 7 + V | V — * ( P ~ ^ — V — + E r + 9 r ■ , D \\\ - . (3.3) 



Moreover, 

ln(n + i) 2 1 
~2 



= - ln(2m) + - In f 1 + — 
2 V ; 2 V 8™ 

1 1 °° I 

= - ln(2m) + - Vf-l)'- 1 — -=-. 

1=1 

Substituting the right-hand side of this last equation into (13. 3p . we obtain 

P — 1 \ 1 



if. = i in( 2m ) + i Et-D'-y^ + 7 + E{E ^(-ir-( 

7=1 p=l k s=0 V 



p — s J 8 P s \ mP 



+ e r + E r + 9 r 



(n + 2 , 



iin(2m )+7+ ^{(-ir^ + x: 

p=l ^ s=0 



^(_l)l>-«/^ _ ^ 1 



2 v ' ^-fl 2n8P ^ 2 s y ' \p-sJ8P~ s \ 

p=l v s=0 v ' J 

+ e r + E r + e r . in °\+l +2 . 

Therefore, we have obtained Ramanujan's expansion into powers of — , and the coefficient 
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Therefore, the formula for H n takes the form 

H n = \ H2m) + 7 + £ { 1 + E (f) M)'^® } -^ + Sr, (3.5) 

p=l ^ s=l ^ ^ ^ 

where 

£ r := e r + £ r + r • ^ + y 2r+2 - (3-6) 

We see t/iat (13. 5p is the Ramanujan expansion with the general formula as given in the 
statement of the theorem, while (13.61) is a form of the error term. 
We will now estimate the error, (13. 6p . 

To do so, we will use the fact that the sum of a convergent alternating series, whose terms 
(taken with positive sign) decrease monotonically to zero, is equal to any partial sum plus a 
positive fraction of the first neglected term (with sign). 

Thus, 
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Moreover 
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where < 5^ < 1 for k = 1,2, ...,r + l and < A r < 1. Thus, £/ie error zs egwa/ to 
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2(r 

= G r • R r +i, 

by fll.6p . where < O r < 1, which is of the required form. This completes the proof. □ 

The origin of Ramanujan's formula is mysterious. Berndt notes that in his remarks. Our 
analysis of it is a posteriori and, although it is full and complete, it does not shed light on 
how Ramanujan came to think of his expansion. It would also be interesting to develop an 
expansion for n\ into powers of m, a new Stirling expansion, as it were. 

References 

[1] M. Abramowitz and I. A. Stegun, Handbook of Mathematical Functions, Dover, New 
York, 1965. 

[2] B. Berndt, Ramanujan's Notebooks, Volume 5, Springer, New York, 1998. 

[3] T. J. l'A. Bromwich, An Introduction to the Theory of Infinite Series, Chelsea, New 
York, 1991. 

[4] E. Cesaro, "Sur la serie harmonique", Nouvelles Annales de Mathematiques (3) 4 (1885), 
295-296. 



13 



[5] Ch.-P. Chen and F. Qi, "The best bounds of the n th harmonic number" , Global J. Math, 
and Math. Sci. 2 (2006), accepted. "The best lower and upper bounds of harmonic 
sequence", RGMIA Research Rep ort Collection 6 (2003 ), no. 2, Article 14. "The best 
bounds of harmonic sequence", |arXiv:math.CA/0 306233, Jiaozuo, Henan, China, 2003. 

[6] D. DeTemple and S.-H. Wang, "Half-integer approximations for the partial sums of the 
harmonic series", J. Math. Anal. Appl., 160 (1991), 149-156. 

[7] J. Edwards, A Treatise on the Integral Calculus, vol. 2, Chelsea, New York, 1955. 

[8] A. Lodge, "An approximate expression for the value of H 1 h • • • H — ," Messenger 

2 3 v 

of Mathematics 30 (1904), 103-107. 
[9] L. Toth, and S. Mare, Problem E 3432, Amer. Math. Monthly, 98 (1991), 264. 
[10] M. Villarino, "Ramanujan's approximation to the n th partial sum of the harmonic se- 



ries" . 



arXiv:math.CA/0402354, San Jose, 2004. 



[11] M. Villarino, "Best bounds for the harmonic numbers", arXiv:math.CA/0510585, San 
Jose, 2005. 



14 



